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Thermoelectric efficiency is defined as the ratio of power delivered to the load of a device to the rate 
of heat flow from the source. Till date, it has been studied in presence of thermodynamic constraints 
set by the Onsager reciprocal relation and the second law of thermodynamics that severely bottleneck 
the thermoelectric efficiency. In this study, we propose a pathway to bypass these constraints using 
a time-dependent control and present a theoretical framework to study dynamic thermoelectric 
transport in the far from equilibrium regime. The presence of a control yields the sought after 
substantial efficiency enhancement and importantly a significant amount of power supplied by the 
control is utilised to convert the wasted-heat energy into useful-electric energy. Our findings are 
robust against nonlinear interactions and suggest that external time-dependent forcing, which can 
be incorporated with existing devices, provides a beneficial scheme to boost thermoelectric efficiency. 


The on-going advances of nano-structure engineering 
have re-energized the search for high-efficiency thermo¬ 
electric devices [H-Q. Till date, almost all studies on 
thermoelectricity are focused on finding high efficiency 
materials guided by the near equilibrium thermodynamic 
quantities like the Seebeck coefficient and the thermo¬ 
electric figure of merit ZT. These passive searches 
have reached saturation and the thermoelectric efficiency 
achieved thus far is still insufficient from a practical 
standpoint^. This is primarily because, in the near equi¬ 
librium regime, the thermoelectric efficiency is limited by 
various thermodynamic constraints, namely, the second 
law of thermodynamics which imposes an unavoidable 
entropy flow and the Onsager reciprocal relation that 
connects the Seebeck and Peltier effects. 

In order to achieve high thermoelectric efficiency an 
active approach to overcome these thermodynamic ob¬ 
stacles is the need of the hour. A possible mechanism 
overcoming these thermodynamic constraints is to apply 
a time-dependent forcing to drive the system far from 
equilibrium. Unlike bulk materials, many nano-systems, 
such as quantum dots d, B , sing le-electron-transistors 
0 , and molecular junctions | 8 l-[l^. can strongly interact 
with an externally applied control force. These systems 
have been the subject of intense theoretical investiga¬ 
tions to better understand the mechanisms underlying 
electron [I^ and heat 0[il transport in presence of 
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a time-dependent control. Several applications such as 
overall device efficiency [iB| . thermopower m , thermal 
refrigeration [T^ , electron pumping [T^ , and heat pump¬ 
ing [i^ [13 have also been studied in these systems to 
figure out the role of an external control. Despite these 
advances the study of dynamic thermoelectric efficiency 
has been highly non-trivial due to the breaking of ther¬ 
modynamic constraints and the non-unique definitions of 
the thermodynamic quantities such as the Seebeck coef¬ 
ficient and the figure of merit ZT. 

In this study we propose a dynamic theory of ther¬ 
moelectric efficiency to overcome the present thermody¬ 
namic limitations such as the Onsager reciprocal rela¬ 
tion. The main idea being that a strong external time- 
dependent control breaks the time-translational invari¬ 
ance of the system and pushes it in the far from equi¬ 
librium regime. This in turn causes a breakdown of the 
celebrated Onsager symmetry relation that allows for the 
possibility to boost the resulting efficiency. The boost 
can be as large as four times the near equilibrium value 
and vitally a large fraction of the supplied input energy 
is constructively utilized to enhance the thermoelectric 
performance. Thus, our novel approach makes available 
an extra knob to engineer high thermoelectric efficiency 
in nano-devices. 

Results 

Dynamic theory of thermoelectricity. Since in the 
far from equilibrium regime the Seebeck coefficient and 
the figure of merit ZT are ill-defined we establish a 
thermoelectric formalism based on the underlying time- 
dependent currents. This objective can be achieved 
through the evaluation of the Onsager transport ma¬ 
trix which relates the electron or heat current with the 
temperature or chemical potential bias. In the conven- 
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tional, near equilibrium, formalism the transport-matrix 
coefficients are autonomous and constrained by various 
thermodynamic laws. However, a time-dependent con¬ 
trol pushes the system far from equilibrium and results 
in the transport coefficients depending on the entire his¬ 
tory of the applied protocol. Moreover the presence 
of a time-dependent control force causes the charging 
and discharging of the nano-system resulting in a time- 
varying thermoelectric current (known as the displace¬ 
ment current) 0123 , which will vanish exactly in the 
near-equilibrium scenario. 

In order to construe the above mechanism we consider 
a two-probe transport set-up consisting of a system con¬ 
nected to a left and right lead as depicted in the inset of 
Fig. [Ha). When the system is subjected to an external 
driving force F{t) the left lead electron and heat current 
will be functions of the thermodynamic forces and the 
entire history of the applied protocol, F(t'), to < t' < t, 
with the starting time tg of the force protocol of otherwise 
arbitrary strength. In the linear response regime for the 
thermodynamic forces, namely, temperature difference 
AT/T and chemical-potential difference Afi/e small, we 
Taylor-expand the currents at each instance of time as, 

Ie(t) = Ieit)\Z=0 + L,,[F]^ + ( 1 ) 

Ikit) = + i2i[F]^ + L,,\F]^. (2) 

Above, the coefficients Lii[F] = Li2[F] = 

d^AT%) Iam=o, L2i[F] = |at=o and L22[F] = 

d(AT/T) indirectly depend on the entire history of 

the applied protocol via the currents 1 ^( 1 ) and I^{t). 
The zeroth order term I^{t) = (OIa^Aq arises only 
in the presence of an external control and is called the 
displacement current. Since the displacement current 
doesn’t change sign under reversal of thermodynamic 
forces we can write the above Taylor expansion as a 
transport-matrix equation that reads, 



(L,,[F] L, 2 [F] Cn-h 

[L2^[F] L22[F] C^[-]) 


( A“^/e\ 
A“T/r 

V nt')) 


( 3 ) 


The super-script a takes values L and R correspond¬ 
ing to the leads. The elementary charge e > 0 and 
A^H = -A^/r = yLh — (similar interpretation for 
A“T). Above represents displacement current 

kernel acting on the history of the applied protocol such 
that the displacement current (t) = [F]. 

For an undriven system; i.e. F{t) =0 Vt, the displace¬ 
ment currents vanish, leaving only the biased currents, 
yielding a nonequilibrium steady state. Specifically, the 
transport matrix reduces to a commonly known, time- 

independent 2 x 2 Onsager matrix L = ( and 

\F21 F 22 ) 

the transport coefficients are obeying the constraints of 


near equilibrium thermodynamic steady state transport; 
namely the Onsager reciprocal relation are valid, impos¬ 
ing that T21/T12 = I. Likewise, the second law of ther¬ 
modynamics ensures a positive thermal conductance, or 
det(L) >00. 

The primary result in this work is to obtain the trans¬ 
port coefficients under a time-dependent control. This 
can be achieved in the following manner: (i) We as¬ 
sume small thermodynamic forces for the temperature 
bias and the potential difference so that the relationship 
w.r.t to these forces stays linear, (ii) The currents are 
evaluated (see below) at any time instant t as a function 
of the two small thermodynamic forces, (iii) Then, set¬ 
ting A“^/e = 0 the slope of the electron- (heat-) current 
w.r.t A“T/T yields Li2[F] [L22[F]) at the time instant t. 
Likewise, for A“T/r = 0 we extract Lii[F] and L2i[F], 
respectively. The intercept of the electron (heat) current 
at time instant t w.r.t A“T/T = 0 or A“/x/e = 0 yields 
the contribution of the currents solely arising from the 
arbitrary driving F(t); i.e. the displacement current. 

In order to investigate the consequences of the time- 
dependent control on the thermoelectric efficiency we 
bias the system with a temperature difference AT, con¬ 
nect a load of resistance i?L to the system and calcu¬ 
late the amount of power consumed by the load. We 
assume that the load is a pure resistor that cannot lead 
to charging effects due to the passage of electron current. 
Therefore, the amount of current passing through the 
load is related to the bias and the transport matrix L. 
After accounting for the back-action from the load, the 
biased electron current reads Ie{t) = Ti2AT/[T(l-f M)], 
where M = Rj^/Rm is the ratio of the resistances with 
Rm = being the resistance of the system. Hence 
the thermoelectric efficiency ratio of the heat-work con¬ 
version reads [2l|, [l^l , 

, . _ _ leRh _ , . 

^ “ det(T)i?M AT/T + L2iRMle - IIRm /2 ' ^ ’ 

Here, we have suppressed the explicit time-dependence 
in all terms on the r.h.s. for notational simplicity. The 
numerator /gi?L is the useful power on the load while the 
denominator is the heat extracted per unit time from the 
hotter lead. The extracted heat consists of three con¬ 
tributions due to the entropy flow det{L)RMAT/T, the 
Peltier heat due to the electron current L2iRMle, and the 
Joule heating term I^Rm with the factor — 1/2 indicat¬ 
ing that half of the heat flows back to the hotter lead. In 
the nonequilibrium steady state this efficiency ratio will 
reduce to the standard formalism 0 where, det{L)RM/T 
represents the thermal conductance, L21RM is the Peltier 
coefficient, and the efficiency is directly related to the 
figure of merit ZT provided that the Onsager reciprocal 
relation L21/L12 = 1 is satisfied. Since all the transport 
matrix coefficients [forming the numerator and denomi¬ 
nator of Eq. (| 4 ])] are affected by the time-dependent con¬ 
trol, a priori it is not clear if the control will have an 
overall enhancing or diminishing effect on the dynamic 
thermoelectric efficiency. Equation (jj]) represents one of 
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our main results, generalizing the conventional thermo¬ 
electric theory to a dynamic one which can be applied to 
far-equilibrium non-steady-state regime. 
Non-interacting electrons. As a proof of concept, we 
first consider a single electron quantum dot in the regime 
of strong Coulomb blockade. The time-dependent exter¬ 
nal force F{t) causes charging and discharging on the 
system and the total Hamiltonian reads, 

H = HL+Hii + Hs{t)+HT, (5) 

where Ha = a = L,R is the Hamiltonian of 

the leads, Ht = J2a=L,RJ2k€a is the tun¬ 

nelling Hamiltonian between the quantum dot and the 
leads, and the Hamiltonian of the quantum-dot system 
is 


Hs{t)=[eo + Fit)]dU. (6) 

This quantum resonant model has been extensively stud- 
ied in the context single-electron-transistors JE S l23l - l^ , 
molecular junctions [2^ or nano-wires [^. The en¬ 
ergy level of the dot can be controlled either via a time- 
dependent gate voltage , or via long-wavelength elec- 
tro mag netic fields such as microwaves [illii or lasers 

The result of heat-work conversion efficiency ratio of 
this model under step-like control is shown in Fig. [1] 
From Fig. [Ua) we detect large enhancements in the ef¬ 
ficiency, upto a factor of 4 compared to the steady-state 
efficiency, as soon as the step-pulse is applied. After some 
relaxation time the values eventually saturate to the new 
steady state. Interestingly, the magnitude of A 21 /T 12 
[Fig. mb)] shows a profile similar to the efficiency ratio 
indicating that the breakdown of the Onsager reciprocal 
relations A 21 /A 12 ^ 1 and the efficiency enhancement are 
closely intertwined. Physically, when L 21 is not bounded 
by L 12 , the contribution of the particle flow to the heat 
current can increase under the influence of external driv¬ 
ing. As a result the efficiency is boosted via increasing the 
useful heat (due to particle flow) while limiting the waste 
heat (due to entropy flow). To substantiate this claim we 
plot det(L) in Fig. (He)- Because the det(T) is propor¬ 
tional to the entropy flow we see that it decreases in the 
regime of efficiency enhancements. Importantly, for suffi¬ 
ciently strong driving (s = O.SF) we detect a regime with 
negative values for det(L), indicating a reversal of the 
entropy flow, even though the overall heat current still 
flows from the hot lead to the cold one. 

Harvested power. Supplementary to the colossal 
boosts in the heat-work conversion efficiency it is also 
crucial that most of the input power due to the control 
force is properly utilized to enhance the heat-work con¬ 
version efficiency. We analyse this using the harvested 
and input power 

h^harvest — (0 (^o) 

winput = -2A(t)/f (t)/e, (7) 
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FIG. 1. Efficiency and transport-matrix coefficients 
for non-interacting electrons, (a): time evolution of the 
thermoelectric efficiency rj (normalized by steady-state effi¬ 
ciency 77 ( 0 ) = 77ss). Inset shows a typical set-np studied in 
this work of a multi-level system (depicted by coloured rings 
within a central sphere), acted upon by an external time- 
dependent control. (b): the entropy flow direction determined 
bydet(L). (c): the Onsager reciprocal relation L 21 /T 12 . The 
control F{t) = s6{t — 1), ksiTL + Tr)I2 = O.IF, chemical 
potential rl = Rr ~ 0, and electron energy eg = 0.5r. The 
efficiency ratio is calculated with a bias ksAT — 0.02r and a 
load resistance 7 ?l = 15hje^. 


where ^harvest is the power harvested from the enhance¬ 
ment of efhciency that is defined as the difference of the 
useful power consumed on the load under time-dependent 
control (first term) and the power from the steady-state 
contribution (second term). Winput is the input power 
via the external control force, where the factor, —F{t)/e, 
represents the input voltage and 21^ {t) is the total dis¬ 
placement current. 

Figure Ha) depicts that the harvested power can be 
much larger than the input power due to the control 
forcing. This feature occurs, even in the linear response 
regime, in the system-parameter regime when the steady- 
state efficiency is low, due to the low electron conduc¬ 
tance, but the Seebeck coefficient itself remains large. 
Thus, the presence of time-dependent control construc¬ 
tively facilitates the movement of electrons and boosts 
the thermoelectric efficiency. 

Resistor-capacitor model. In order to better un¬ 
derstand the displacement current in the high temper¬ 
ature and weak control limit, we propose an elemen¬ 
tary resistor-capacitor model as illustrated in the inset 
of Fig.Hb)- The time-dependent control is acting on the 
gate with capacitance Cg which can induce charging or 
discharging of the capacitor. This variation leads to a 
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FIG. 2 . Power harvested, displacement current, and 
the intuitive resistance-capacitance model, (a) the 

time-dependent harvested power (red dashed line) and the 
input power due to driving (black solid line). The system is 
a non-interacting electron model with load Rl = 50ft/e^ and 
biased with ksiTL -|- Tii)l2 = IF, AT = 0.2r, /j-l = fJ-R = 0. 
(b) displacement current (solid line) and the fit using the re¬ 
sistance capacitance model (red circles) for non-interacting 
electron model with Tl ~ Tr = IF, fiL = fJ-R = 0, and 
to = Ift/r. The fitting parameters are R = 11.6fi/e^ and r = 
1.06/i/r. The common system parameters are: s = O.OOIF 
and £o = 2.5r. 

current generation which flows from the capacitance to¬ 
wards the leads which are represented as the two sink 
sources (the ground connection) in the circuit. The cur¬ 
rent generated solely depends on the time-dependent con¬ 
trol and does not require a thermodynamic bias between 
the leads for its existence and is known as the displace¬ 
ment current (t). Due to dissipative effects the current 
experiences a total resistance R while flowing from the 
capacitor to the leads. 

In case of the quantum dot model with non-interacting 
electrons subjected to a step-like gate control F{t) = 
s6{t — to) the solution reads, 

= ( 8 ) 

The intuitive picture for the displacement current above 
is based solely on circuit law considerations. Hence a 
priori it is not clear if such a model is able to describe 
correctly a fully quantum mechanical system. In order to 
justify that this indeed is the case we use the parameters 
R and r from Eq. ([8]) as variables and fit the equation to 
the fully quantum mechanical displacement current ob¬ 
tain via nonequilibrium Green’s function. Figure [D(b) 
shows the NEGF calculation as a solid line and the fit 
via the red dots. The perfect fit gives us the param¬ 
eters T Ri h/T which further strengthens our resistance- 
capacitance circuit model. This is because in an open dis¬ 
sipative quantum system one expects the relaxation time 
of the system to be inversely proportional to the sum of 
the coupling strengths of each lead T Thus, the 

verification of our resistance-capacitance model gives an 
intuitive picture of the displacement current. 
Electron-phonon interaction. One of the main chal¬ 
lenges for experimental devices to obtain enhanced effi¬ 
ciencies is the unavoidable presence of nonlinear inter¬ 


actions mainly arising due to the involvement of stray 
phonon modes [s^, which can drastically alter the ther¬ 
moelectric efficiency. Hence, we consider a quantum dot 
interacting with a single phonon mode giving rise to the 
following electron-phonon interaction Hamiltonian, 

Hs{t) = [eo + F{t)]dU + ojoo^ a + Xd^d{a^ + a) . (9) 

Here, and a are creation and annihilation operators 
of the phonon, wq is the phonon angular frequency, A is 
the electron-phonon interaction strength and Fit) repre¬ 
sents the time-dependent control of the coherently driven 
quantum dot. The model manifests itself under various 
phwical scenarios like in a nano-mechanical resonator 
[ 1 , 13 , [ 13 , molecular junction [H, and standard lat¬ 
tice vibration model [35l| . Recently it was shown that 
a small amount of nonlinearity in this model can greatly 
suppress the steady-state efhciency [s^. Thus, the model 
serves as a perfect test bed to establish the robustness of 
our approach. 

In Fig. [3 we depict the results for the interacting 
electron model. In case of the delta shape and square 
wave driving we modulate the system for sometime and 
then let it relax to reach its nonequilibrium steady state. 
Clearly the enhancement in the efficiency (as seen from 
the bottom row of Fig. El) is observed even for a rela¬ 
tive strong nonlinear interaction A as long as the sys¬ 
tem dynamics is time-dependent. The long-time limit, 
when the transient effects are wiped out, is easily recov¬ 
ered for the delta shaped and square wave driving giv¬ 
ing ri{oo)/r]ss = 1 [not depicted in Fig. E]- Analogous 
to the non-interacting electron model, the enhancements 
are closely related to the breakdown of the Onsager re¬ 
ciprocal relation L21/L12 and the second law of thermo¬ 
dynamics det(T). Interestingly, external forcing alone is 
not sufficient to enhance the systems efficiency as seen 
from the case with a periodic sinusoidal driving where 
the efficiency even decreases when T21/T12 < 1 . Thus, 
we speculate that although driving is a necessary condi¬ 
tion to allow the breakdown of stringent constraints it 
does not sufhce to enhance the efhciency of the device. 
One possible sufficient condition for an enhancement is 
the abrupt variation in the driving held which causes a 
sudden change of the charge occupation in the system. 
As a result the displacement current will be large (fourth 
row of Fig. El). 

Discussion 

The performance of modern thermoelectric devices is 
measured in terms of its efhciency to convert waste heat 
energy into useful electric energy. Till date, most efforts 
towards efhciency enhancement are based on the search 
for suitable materials guided by the near-equilibrium 
thermodynamic constraints. These efforts have led to 
considerable improvements in the heat-work conversion 
efhciency but the quest for commercially feasible efh¬ 
ciency has been futile so far. 

In the present study, we provide an active and compli¬ 
mentary approach to the material search avenue. This is 
achieved by introducing a time-dependent control which 
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FIG. 3. Efficiency and transport-matrix coefficients 
for electron-phonon interaction. Plot of the entropy flow 
represented by det(L) (a, b, and c), the Onsager recipro¬ 
cal relation L21/L12 (d, e, and f), the displacement current 
(g, h, and i), and the efficiency ratio normalized by the 
steady-state r]{t)/ris3 (j, k, and 1) for the interacting elec¬ 
tron model. The system is subjected to delta pulse driving 
-P'(i) = s J2n = {1, 5, 7, 9,11} (a, d, g, 

and j), multi-step driving = s when fit € [1, 3] U [5, 7] 

and F{t) — 0 elsewhere (b, e, h, and k) and a periodic si¬ 
nusoidal drive F{t) = 2s9{t — to) sin(fl7rt) (c, f, i, and 1), 
where fl = lOF/fi and to = G.lh/T. Other parameters 
are ksiTL + Tii)l2 = IF, AT = 0.2r, = fJ-R = 0, 

TL = Fi^ = F/2, £0 = 2F, ujQ = lOF, A = 3F, and s = IF. 


pushes the system far from equilibrium and provides 
a rich playground without thermodynamic limitations. 
The control force can be fused with existing high effi¬ 
ciency devices and would allow us to further boost their 
efficiency by a factor of 4. The enhancements are robust 
and persist even in presence of nonlinear interactions in¬ 
dicating its usefulness in existing experimental set-ups. 

Thus, our work opens up a whole new arena where we 
shift attention away from a material design perspective 
and focus on the non-trivial far from equilibrium regime 
which leads to smart device design. Overall the method 
presented herein provides a rigorous stepping stone which 
can have wide ran ging im pact for fields such as ther¬ 
moelectric cooling [36l - l39l |. solar thermoelectric genera¬ 
tion (inl - E^ . and can be extended to the field of spin 
caloritronics to efficiently pump spin-currents using ther¬ 
mal gradients (d^ . 

Methods 

Nonequilibrium Green’s function. For non¬ 


interacting electron in a quantum dot under step-like 
control with F{t) = s0{t — to) within the wide-band ap¬ 
proximation, i.e., ra(e) = |V)“p(5(e: - Sk) = r/2 

(a = L,R), an exact solution of the electron and heat 
currents can be obtained using the Landauer formalism 
via the NEGF approach , reading 

/ oo J 

-/„,(e)iF,7,)(e,t), (10) 

a=n,n. 

where the kernels 


Kr'{e,t) = eZ^^'{e,t), 

Kr' ie,t) = ^lm{A{e,t)dtA*{e,t)} 
-b (e - 


Z° 




6a,aAm{A{e,t)} + -\A{e,t)\'^ 


( 11 ) 


The Fermi-Dirac distribution of the a-th lead /a(e) = 
[1 -(-with /3a = l/iksTa) and A{e,t) is the 
spectral density, 

“ (£ - £0 + *r/2) {e-eo-s + tT/2) ^ ’ 

Using Eq. (I10|l it is possible to obtain the displacement 
current. 


^ / de^lm[A{e,t)] + ^\A{e,t)\^]f{e), 

(13) 

by setting /L(e) = fnis) = /(e)- For step-like driving 
we obtain the intuitive results where the displacement 
current vanishes in the long time limit. Subtracting the 
displacement currents from the total currents we obtain 
the biased currents. 


r2 r°° rjp 

Ie{t) = ^J ^^\A{e,t)f[fL{s)-fR{e)], (14) 

■p pCiO J 


Im 


A/e, t)h 


dA*{e,t) 
~dt 


}[/L(e) - /fl(e)].(15) 


The biased electron current I^/t) and the first term of 
the biased heat current /jf (t) take the form of the Lan¬ 
dauer formula which helps preserve the Onsager sym¬ 
metry. The second contribution to the biased heat cur¬ 
rent arises only due to the presence of an explicit time- 
dependent, time-reversal breaking control that is respon¬ 
sible the breakdown of the Onsager symmetry. 

Given the biased and displacement currents the NEGE 
formalism allows us to obtain a closed form expression 
of the efficiency rj{t) in the weak system-bath coupling 
and weak control limit. In order to achieve this goal we 
first simplify the efficiency given by Eq. (|3|) for maximum 
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power output, i.e., M = 1 and negligible Joule heating 
(small AT) to obtain 


v{t) 


Lh at 

4Liii22 ~ 2L12L2I T 


(16) 


The weak coupling limit (see Append, of Ref. [i^) trans¬ 
forms the transmission |A(e,t)p to a delta function and 
subsequently keeping only leading order terms in the con¬ 
trol strength s helps obtain the currents analytically in a 
closed form. These currents are then used to obtain the 
transport matrix coefficients Lii,Li2,T2i, and L22 that 
result in the efficiency 


Vii) ='nss + s 


dVss f2{t - to) 
deo fi{t-to) 


, at 

-b sr]{t - to) — , 


(17) 


The protocol F{t) begins at to (0 < to < t) and ends 
at time t and the displacement current depends on the 
complete history of the protocol. 

Quantum master equation. Due to the presence of 
nonlinear interactions in the system we resort to the time- 
dependent quantum master equation formulation to eval¬ 
uate the currents. The formulation treats the nonlinear 
interactions exactly under an arbitrary forcing at the cost 
of a weak system-lead coupling. Following the standard 
scheme [l^ the quantum master equation for the reduced 
density matrix p{t) of the system reads 

^ 1: RSh. p«. (22) 

k,k' 


where rjss is the steady-state efficiency at time t = 0 and 
flit — to) can be expressed in the high temperature limit 
as, 


Vit - to) 


TJ" {t — to)e 
2h{p,-eo) /i(t-to) 


(18) 


The two functions are, /i(t — to) = -b /2(t — to) 

and /2(t — to) = (1 — Their ratio, /2(t — 

to)/fi{t — to), monotonically increases from 0 to 1 as 
(t — to) —>■ 00, thus providing no temporal boost in the 
efficiency arising from the second term in Eq. (flTl) . The 
third contribution, i.e., oc fi{t — to), arises due to the con¬ 
tribution that breaks the Onsager symmetry in the biased 
heat current [second term of Eq. (ITS]) ] and is responsible 
for the temporal boost in the thermoelectric efficiency 
(in the leading order of s). This further strengthens our 
claim that it is indeed the breaking of Onsager symmetry 
leads to a boost in thermoelectric efficiency. 

Equations of motion for the resistor-capacitance 
model. We further elucidate on the resistor-capacitor 
model as shown in the inset of Fig. EJb). Consider that 
the capacitor has a charge Q then the voltage on its upper 
plate will be the sum of the voltage across the resistances 
and the voltage across the capacitor [d^, namely. 


m 

e 


Q{t) 

Co 


2/f(t)i?. 


(19) 


Above since R is the total resistance, 2R will be the re¬ 
sistance of each resistor giving the voltage across each 
resistor as 2I^{t)R. Differentiating the above equation 
with respect to time we obtain 

i^[t) + h^{t) + ^F{t) = t), ( 20 ) 

r 2eR 

where r = RCg represents the relaxation time of the 
leads. Above since the displacement current is due to 
the charging or discharging of the gate capacitance Cg 
we have used Q[t) = 21^ {t) as the total displacement 
current. The solution to the differential equation reads 


/f(t)=£f[T] 


1 


dt'F{t')eC-t)/r^ 


( 21 ) 


where the relaxation four-tensor 


7^: 


ijk 


(t) - YnlYlf (^) 


-b C.C.. 


Above Aij{t) = Ei{t) — Ej{t) is the energy spac¬ 
ing with Ei{t) as the i-th instantaneous eigenenergy. 
Since the external time-dependent driving only modu¬ 
lates the eigenenergies, and does not affect the eigen¬ 
states of the system Hamiltonian, we use the eigenstates 
of the static Hamiltonian as our basis. The transition 
matrix Wif'(<) = e-^t' /h^jkk' 

with the correlation function (t) = {B^{t)B^ (0)). 
The vector-operators Y and B belong to the sys¬ 
tem and lead Hilbert space and appear in the 
tunnelling Hamiltonian; i.e., Y = {d, rff} and 

B = {Z)a=L.fl Sfcea S)a=L.fl Sfcea 

yfe (^B^) denoting the fc-th component of the Y (B) vec¬ 
tor. The operator B{t) is the free-evolution of B with 
the lead Hamiltonian + i7_R. 

Generalizing the nonequilibrium steady-state formula¬ 
tion [H, m, to encompass time-dependent control 
E(t) we obtain the expression for currents as. 


^ E (<)] } , (23) 

k,k' 


where the electron or phonon hopping rates 
are defined, similar to the master equa¬ 
tion, using the current-lead correlation func¬ 
tions C^p[)(t) = (H''(t)Hg(^)(0)), where the op¬ 

erators Be = {eX:feeL^fc^4>-eEfe^fc^Cfe} and 

-ML)R/'Cfc }. 

Once we know the time-dependent reduced density 
matrix of the system, we can deduce the displacement 
current as the time-derivative of the average charge on 
the dot as, leit) = —Q/2, where Q = Tr[p(t)dfd]. 

The above quantum master equation formalism is valid 
in the weak system-bath coupling limit and holds true for 
arbitrary control strength s and electron-phonon interac¬ 
tion strength A. Moreover, we do not resort to the secular 
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(or rotating wave) approximation and the master equa¬ 
tion is kept seini-non-Markovian since time t is explicitly 
present in the integral-limits of the transition matrix. 


The approach is robust to deal with nonlinear interac¬ 
tions exactly and thus allows to obtain the transport- 
matrix coefficients for strongly nonlinear systems. 
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